Quantum simulation of small-polaron formation with trapped ions 
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We propose a method for simulating polaron physics using a one-dimensional system of trapped 
ions acted upon by off-resonant standing waves. This system, envisioned as an array of ion micro- 
traps, in the single-excitation case provides a realization of the anti-adiabatic regime of the Holstein 
model. We show that the strong excitation-phonon coupling regime, characterized by the formation 
of small polarons, can be reached using realistic values of the relevant system parameters. Finally, we 
propose measurements of the quasiparticle residue and the average number of phonons in the ground 
state, experimental probes validating the polaronic character of the phonon-dressed excitation. 

PACS numbers: 03.67.Ac, 37.10.Ty, 71.38.Ht 



The field of quantum simulation inspired by the 
ideas of Feynman Q , holds promise to advance our un- 
derstanding of complex many-body systems Q. In this 
context, trapped ions Q constitute a versatile experimen- 
tal platform on which to explore a variety of models and 
phases [f|. Owing to the high degree of control over the 
internal degrees of freedom and the single-ion address- 
ability, trapped-ion systems allow both analog and digi- 
tal quantum simulations of quantum spin models [6rl8], 
models with bosonic degrees of freedom [9|, and even of 
phenomena outside the realm of low-energy physics [Toj j . 

Electron-phonon (more generally, particle-phonon) 
coupling, a traditional subject of solid-state physics, has 
only quite recently attracted interest from the quantum- 
simulation community The most striking conse- 



quence of this coupling in semiconductors and insula- 
tors is small-polaron formation: an excess electron (or a 
hole) can be severely localized in the potential well that 
it creates by displacing the surrounding atoms of the host 
crystal ( "self-trapping" ) [l2| ■ While the polaron concept 
was conceived by Landau and Pekar in their studies of 
polar semiconductors evidence has by now accumu- 
lated for small polarons in systems as diverse as amor- 
phous and organic semiconductors, manganites, undoped 
cuprates, superlattices of graphene, etc. [14j |. In particu- 



lar, studies of small polarons in systems with short-range 
electron-phonon coupling are typically based on the Hol- 
stein molecular-crystal model [l5j. The static and dy- 
namical properties of this model, describing purely local 
interaction of tightly-bound electrons with dispersionless 
(Einstein) phonons, have been extensively studied 12[. 

In this Letter, we present a scheme for simulating 
small-polaron formation using a linear array of trapped 
ions subject to off-resonant standing light waves. We 
show that, if the relevant parameters are chosen appro- 
priately, this system is described by an effective model 
whose single-excitation sector corresponds to the anti- 
adiabatic regime of the Holstein model (hopping ampli- 
tude much smaller than the phonon frequency). Fur- 
thermore, we find that with realistic values of experi- 



mental parameters (laser intensities, detunings, trap fre- 
quencies, etc.) a regime is realized where the dressed 
single-excitation ground state has polaronic character. 
Finally, we suggest methods for measuring the quasipar- 
ticle spectral residue and the average number of phonons 
in the single-excitation ground state, quantities that pro- 
vide typical signatures of small-polaron formation. 

System and effective Hamiltonian - We consider a one- 
dimensional system of N trapped ions with four internal 
states (|t) z > l e t) Z ' \ e i) z ) an d a possible physical 

implementation as an array of microtraps (for an illus- 
tration, see Fig. [J). We take the z-axis to be parallel 
to the direction of the ion chain (longitudinal direction) , 
while x and v will be the radial (transverse) directions. 
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FIG. 1: (Color online) System schematic: Trapped-ion chain 
subject to three pairs of counterpropagating laser beams giv- 
ing rise to three standing waves. The standing wave with 
the polarization a + in the direction a induces the transi- 
tion from |f) Q to \ei) a , where \t,i) x = (|t>* ± I l) z )/V2 and 
tt,l) y = (tt) z ±i\l) z )/V2. 
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The total potential felt by an ion consists of the 
trapping-potential contribution (with trapping frequen- 
cies uj a , a — x,y,z) and the Coulomb repulsion between 
the ions. In the harmonic approximation, the total vibra- 
tional Hamiltonian of the system reads (hereafter h = 1) 

The displacements qf of the ?-th ion (mass m) from its 
equilibrium position zf (pf are the corresponding mo- 
menta) are expressed through the creation and annihi- 
lation operators of phonon modes with frequencies u SiCt 
as it = E s - A/f kQ( 2TOW S:a)~ 1/2 ( a ^a + a\ a ), such that 
-frph = Yl s ,a u s ,eta\ a a Sta ; M l s a is the z-th component 
of the normalized eigenstatc M. s a of the matrix JC a : 

Yji,j - v l- n^ 7 ; - v ^ „ = Us,a s aa>- If the trapping potential 
in a transverse direction a dominates over the Coulomb 
interaction [the stiff limit: j3 a = e 2 / (rruJ^d^) -C 1, where 
do is the mean distance between two adjacent ions], the 
vibrations in this direction are highly localized in real 
space (weakly dispersive phonons) p. On the other 
hand, by tuning uj z one can reach both the weakly- 
dispersive (fi z <C 1) and the strongly-dispersive [fi z 1) 
regime for longitudinal phonons. 

As shown in Fig. [1] three pairs of counterpropagating 
laser beams with the polarization cr + form three standing 
waves, which induces the transition from |f) to |e^} [f|. 
Under the rotating-wave approximation, the adiabatic 
elimination of the higher energy states |e-f) and \e±) 
leads to the effective interaction 

^ = E^ cos2 m + ^)( 1+ <) ( 2 ) 

2, a 

between phonons and the pseudo-spins corresponding to 
the internal states |f) z and \l) z , where A a = Wi jQ — ojq 
is the detuning between the energy-level spacing ujq and 
the laser frequency along the a direction; of is the Pauli 
matrix in the basis of \\) z and \\.) z . These pseudo-spins 
are simultaneously acted upon by global magnetic fields, 
as described by the Hamiltonian H m = ^\ B a erf . 

In the Lamb-Dicke limit, the ion-position dependent 
coupling constants can be linearized around the equi- 
librium position, so that the interaction term becomes 
H i = Ei.a^gf (1 + erf). Here, the effective force F a 
is given by F a ~ G^/c Q /(2A Q ). For convenience, we fo- 
cus on the parameters F x — F y = F, w x = uj y = ujq, and 
Px = Py = P- Given that Fj {<jj ^J2muj ) = n <C 1, due to 
the large transverse-phonon frequencies, one can pertur- 
batively eliminate transverse phonons using a Frohlich- 
like canonical transformation [6]. In this manner, we ob- 



tain the effective Hamiltonian 

H eS = ff L +i^ 9 ?(l + o?) (3) 

i 

+b z 2^°i+ 2^ , j » 

i i<j,a—x,y I J I 

where the external magnetic fields B a = F 2 /muj 2 t (a = 
x, y) were chosen in order to cancel the effective fields 
from the adiabatic elimination. The Hamiltonian 

^ = - 2 E^-Err^(^+ H - c -) ( 4 ) 

describes longitudinal phonons in terms of local-phonon 
operators, where qf — {2mio z )^ 1 / 2 {b i + b\) and the 
renormalized longitudinal-phonon frequency is given by 

£ 2 = w^l + Ej'^A/l*-/! 3 )- The form of the sec- 
ond term in Eq. results from omitting the phonon- 
number nonconserving terms under the rotating-wave ap- 
proximation; the corresponding condition coincides with 
the stiff limit j3 z = e 2 /(mw 2 z dl) < 1. 

By switching to the fermion representation of the ef- 
fective spin operators [I + erf — > 2c\ci, crfcr* + v\o V j — > 

2(c\cj + He.)], the effective Hamiltonian in the single- 
excitation sector of the problem can be written as H se = 
H c + H L + He-ph, with 

H ° = E 7-^3^ + H - c -) + 2B * E c ^ > 

He- P h = ffcZ) z y^c|ci(bj + &|) . (5) 

i 

Here J = A/3t] 2 ujq is the hopping amplitude, while 
g = \[2F z j \JmCb\ is the dimensionless excitation-phonon 
(e-ph) coupling strength. This coupling, described by 
-ffc-ph, has Holstein-like (local) form. 

Limits of validity and parameter range - Let us sum- 
marize the conditions of validity of the Hamiltonian H sc 
and specify in which parameter regime it reduces to the 
Holstcin model with dispersionless phonons. 

Firstly, the use of perturbation theory for eliminating 
transverse phonons requires that r\ <C 1, which implies 
that J -C wo must be fulfilled. Secondly, while a weak 
phonon dispersion [here proportional to f3 z u} z ; cf. Eq. (j4|] 
is inherent to the stiff limit, for our purposes the disper- 
sion should be negligible even compared to the hopping 
amplitude J. From the condition f3 z Co z <C J we readily 
obtain loq/ui z < 4?y 2 < I. Thus our system can simulate 
the Holstcin model with J <C wo < &z, i-e., in its anti- 
adiabatic regime. In other words, the effective Hamilto- 
nian of the system then reads H = H c + H p h + iJ e _ p h, 
with H ph = u) z J2i b\bi. 

Unlike the genuine Holstein model, which only involves 
nearest-neighbor hopping, the last effective Hamiltonian 



also contains long-range hopping terms decaying with the 
characteristic inverse third power of the distance. How- 
ever, in the anti-adiabatic regime that we will be solely 
concerned with these long-range terms are of minor im- 
portance. 

We now discuss whether one can realize the strong- 
coupling regime of the Holstein model, characterized by 
small-polaron formation. To achieve that, a coupling 
constant g of the order of unity is needed. We first ex- 
press g in terms of the relevant experimental parameters 
starting from the expression F z ~ (G z /A z ) x (2n/X), 
where A the laser wavelength, and A z the detuning be- 
tween the laser and the internal energy-level spacing, 
bmce guj z = y/2F z /s/muT z and a/d ~ 1/(AV2 muj z ), 
where a is the size of the trap ground state (for con- 
creteness, we assumed that do = A, while in general do is 
an integer multiple of A), we readily obtain 
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FIG. 2: Calculated quasiparticle residue at k = for different 
values of the renormalized longitudinal-phonon frequency ui z ■ 
All three curves correspond to the hopping amplitude J = 
25 KHz. 
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Thus the value of g can be varied by tuning the ratio 
G z / A z , for different values of the phonon frequency uj z . 

We choose A = do = 2fim, and the typical value of 
the detuning A = 1000 GHz. With Rabi frequencies 
G z = 10 - 100 GHz, the ratio G 2 Z /A Z is in the range 
0.1 — 10 GHz. At the same time, the phonon frequency 
is typically u! z /2n = 1 — 20 MHz. For example, already 
for G 2 Z /A Z = 2 GHz and cD z /2tt = 10 MHz (the ratio 
a/do ~ 4.25 x 10 -3 ) we obtain g 1.72, a value of g that 
belongs to the strong-coupling regime. 

Polaron ground state — To determine the ground state 
of the system, we make use of the variational method. 
Variational approaches, commonly employed in small- 
polaron studies, have been established to yield results 
that agree very well with exact-diagonalization results 
on systems of small size 18| . 



The eigenstates of the total Hamiltonian are good 
quasi- momentum states, i.e., simultaneous eigenstates of 
the total momentum operator K — J2k^ c k c k~^J2 q Qb q b q 
(with eigenvalues denoted by k in what follows). There- 
fore, variational states must be Bloch-like linear com- 
binations \tp K ) = N~ x / 2 e tKn \^ K (n)) of Wannier-like 
functions \ip K (n)}. In particular, for the Toyozawa Ansatz 
states 19] 



JV/2-1 

\*l>*(n)) = /2 $ «( m ) e 

m=~N/2 



mm t 



J0> e ®|£ re (n)) ph , (7) 



with |£ K (ra)) ph = IL ex P { v i b i+i - wf*6 n+ i)|0) ph being a 
direct product of phonon coherent states at sites n+l (I — 
-N/2, . . . , N/2). The parameters {$«(m), uf }, in total 
27V of them, are complex valued. The normalization is 
given by (if> K \i> K ) = T, m .m> K(™.')$ K (m) 5^,_ m , where 



S«_ n , =exp[-|Ep(l«p- 



w p—n I 



-2v K 

p—n 



'Vp-n)]- 



In the limit Cb z / J — » oo the Toyozawa Ansatz yields 
(for k — 0) the exact ground state of the Holstein 
model, given by aj l0 |0} e <g> \v*=° = g) ph (where wf =0 at 
sites I ^ no are equal to zero) [l7j]. The last statement 
holds true not only for the genuine Holstein model with 
nearest-neighbor hopping J, but also for our model with 
long-range hoppings, as the latter have amplitudes pro- 
portional to J. This makes the Toyozawa Ansatz partic- 
ularly appropriate for the problem at hand. 

The ground-state energy Eqs is obtained by minimiz- 
ing the expectation value (ip K= o\H\ipK=o) with respect to 
AN variational parameters, the modules and phases of 
the complex parameters <E> K= o( m ) an d vf =0 ' . While there 
is no principal obstacle to treating even larger system, in 
the following we present results obtained for N = 32. 

Characterization of the polaron crossover- The quasi- 
particle residue at quasi-momentum k is defined as 
Z k = \(^ k \ip k }\ 2 /{i>k\ipk), the overlap squared of 
the bare-excitation Bloch state \^k) = c jj0) = 
N^ 1 / 2 J2 n e lfc ™cJi|0) and the corresponding (dressed- 
excitation) Bloch state |V>fc) of the coupled excitation- 
phonon system. For the Toyozawa Ansatz state, we find 



J2 *i(m')**H S^ f . 



-\<\ 2 



(8) 



■m I 



The quasiparticle residue Z k= o, when evaluated for op- 
timal values of the variational parameters, characterizes 
the ground state of the system. The change of this quan- 
tity from unity (non-interacting system) to values very 
close to zero (Fig. [2]) illustrates the crossover from a 
quasi-free excitation to a small polaron. 

Another quantity characterizing the polaron crossover 
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FIG. 3: Average number of phonons in the single-excitation 
ground state. All three curves correspond to J = 25 KHz. 



is the average number of phonons in the ground state 



N ph EE (^ =Q | Ys b i bi l^=°) 



£ ^(m^W^,^ . (9) 



As illustrated in Fig. [3l this quantity varies from 7V p h ~ 
(quasi-free excitation) to values N p h > 3 (small polarons) 
with increasing coupling strength. 

Our analysis indeed shows that the effect of the fast- 
decaying long-range hopping on the quantities charac- 
terizing the polaron crossover is very small in the anti- 
adiabatic regime. This is consistent with a recent study 
of the Holstein model with next-nearest-neighbor hop- 
ping which established that such hop ping is significant 
only in the opposite, adiabatic regime 20 1. 

Experimental validation - We propose two experimen- 
tal methods, as illustrated in Fig. [TJ to detect the quasi- 
particle residue Zk=o and the average number of phonons 
iVph in the single-excitation ground state. 

To probe Z^—q, we first prepare the system in the 
ground state (without a polaronic excitation), namely, 
the vacuum state of the longitudinal phonons and the 
internal states \\) of all the ions. Then, all the laser 
beams are turned on. At the same time, an additional 
Raman laser with the detuning 5 is used for generat- 
ing an excitation (polaron). The latter is described by 
the Hamiltonian ff dc tcct = £iV~ 1/2 EiOi" e ~ i5 *+ H - c -)' 
where £ is the effective Rabi frequency. After the evo- 
lution time T, we switch off the laser beams and measure 
O = ^~ 1 ^ 2 J2i a i m- By repeating the measurement 
one can determine the mean- value (O). For short times 
T -C l/£, near the resonance frequency S ~ Eqs of the 
ground state in the single-polaron space, first-order per- 
turbation theory yields 



(O) ~ 2£Z fc= 



siu(iT) - sin(£; G sT) 
' S-E GS 



Thus measurements of (O) provide a direct access to the 
quasiparticle residue Zk=o, the quantity whose strongly 
reduced values indicate the presence of small polarons. 

For detecting the mean phonon number N p h, we first 
prepare the system in the single-polaron ground state 
|£7gs) = \ipK=o) an d switch off all laser beams. Then, 
we detect the internal state of the ion on the site j. If 
the internal state is •, the Raman laser is turned on 
to induce the transition between to another state 
\si)j] if the internal state is the Raman laser is 

turned on to induce the transition between |f) ■ to the 
state \sf)p as shown in Fig. [TJ Here, each Raman 
laser is turned to the lower sideband detuning S = io z . 
Thus, in the interaction picture the time evolution of 
the system is described by the Hamiltonian H a= ^^ = 



(a\ +H.c). After the evolution time T, we de- 



tect the probability p a j of the transition to \s a )j- For 
the short-time evolution, first-order perturbation theory 
leads to p aJ = i 2 T 2 (E GS \b]b 3 P^ ] \E GS ), where Pj j) = 
\a)j (a\ is the projection operator. Finally, the average 
number of phonons at site j is Nj — YlcrP<r,j/(£ 2 T 2 )- 
Therefore, measurements of the transition probabilities 
p a . j (j = 1, . . . , N) allow us to extract the average num- 
ber of phonons at each site (microtrap) . Large values of 
iVph = Y^j Nj represent a signature of small polarons. 

Conclusions- To conclude, we have suggested a 
scheme for the quantum simulation of polaron physics 
with trapped ions. We have shown that with realistic 
values of the relevant experimental parameters one can 
reach the strong-coupling regime of the Holstein model, 
characterized by small-polaron formation. An experi- 
mental realization of our suggested scheme in an array 
of ion microtraps is clearly called for. The scheme pro- 
posed here opens up new possibilities to analyze polaron 
physics with trapped ions. For instance, one may go 
beyond the single-excitation regime, and thus study the 
interactions among polaronic excitations. Furthermore, 
with the advent of surface traps 0, [2l[ , it should be pos- 
sible to extend the study to two dimensions, and also to 
engineer other kinds of interactions. Experiments with 
trapped ions may allow us to better understand polaron 
physics in these circumstances, and to develop new the- 
oretical tools to describe the resulting phenomena. 
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